Light scattering from disordered overlayers of metallic nanoparticles 
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We develop a theory for light scattering from a disordered layer of metal nanoparticles resting 
on a sample. Averaging over different disorder realizations is done by a coherent potential approx- 
imation. The calculational scheme takes into account effects of retardation, multipole excitations, 
and interactions with the sample. We apply the theory to a system similar to the one studied ex- 
perimentally by Stuart and Hall [Phys. Rev. Lett. 80, 5663 (1998)] who used a layered Si/Si02/Si 
sample. The calculated results agree rather well with the experimental ones. In particular we find 
conspicuous maxima in the scattering intensity at long wavelengths (much longer than those corre- 
sponding to plasmon resonances in the particles). We show that these maxima have their origin in 
interference phenomena in the layered sample. 

PACS numbers; 78.67.-n, 42.82.Et, 42.25.Dd 
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I. INTRODUCTION 

During the last several years there has been an in- 
creased interest in near-field optics. This development 
has been driven in part by the development of the near- 
field scanning optical microscope, but also by a growing 
awareness of the importance of near-field effects and the 
influence they have on the radiated, far fields. Examples 
of probes that investigate or make use of the coupling 
between near-field and far-field effects include surface- 
enhanced Raman scattering (SE]^),El light emission from 
scanning tunneling microscopesp transipission through, 
and reflection from,|-parrow gratings Jju and scattering 
from particle arrays JjQ Applications range from single- 
molecule spectroscopy using SERS to creating microcav- 
ity light sourceaS and using the plasmon resonances of 
metallic nanoparticles for the transmission of "optical" 
signals at a sub-wavelength length scaleB 

The optical properties of small metallic particles and 
collections thepeQf-|tja|We been studied inteiT^vely both 
experimcntaUyli3li3ljT3 and theoretically.E5t3 Many of 
these studies have focused on measuring and calculating 
properties such as average transmission and absorption 
in a colloidal solutions or surface fllms, and to determine 
the role of inter-particle interactions on these properties. 
At present there is also a lot of interest in studying the 
interaction between surface plasmons and particles rest- 
ing on the surface in question since the presenne,p£ the 
particles can induce surface-plasmon band gaps rWA 

In this paper we will focus on, and develop a the- 
ory for diffuse light scattering from disordered overlay- 
ers of nanoparticles resting on a sample. The motivation 
for this work was initially pcovided by an experiment in 
which Stuart and Hall (SII)eI measured the diffuse light 
scattering from random silver nanoparticle arrays fabri- 
cated on top of dielectric as well as semiconducting and 
metallic samples. The silver particles had a a size of some 
100 nm and the intensity of scattered light as a function 
of photon energy showed rather intriguing features. The 



most striking result was that the scattering intensity had 
pronounced maxima at very long wavelengths (A > 1000 
nm) when the sample consisted oLSi and Si02 in a mul- 
tilayer structure. Stuart and HallO attributed this to an 
enhancement of the dipole-dipole interaction between the 
different particles mediated by waveguide modes that are 
supported by the layered sample. In this scenario, inci- 
dent light is initially scattered into a waveguide mode in 
the sample by the nanoparticle array. The electromag- 
netic field associated with a waveguide mode is in general 
evanescent outside the sample. Therefore light scattered 
into such a mode is "trapped" by the sample and can 
travel a relatively long distance before it is damped due 
to losses in the sample or because it is rescattered into 
free space by the nanoparticles. Thus, it is clear that 
the presence of a sample can lead to an increase in the 
particle-particle interactions. Nevertheless, as we will see 
below, the results of this calculation do not support the 
scenario of Stuart and Hall. 

A rigorous theory that can be used to calculate the 
scattering intensity from the system in question must in- 
clude a number of ingredients: particle-particle interac- 
tions through near-field and far-field interactions as well 
as interactions mediated by the sample have to be taken 
into account. In addition, the nanoparticle layer is disor- 
dered and the calculation must deal with the averaging 
over different disorder realizations. Here we do this by 
means of a lattice-gas modeJla in which the sites of a 
regular lattice are either occupied by spheres of a cer- 
tain size and composition or else empty; the averaging 
over different realizations of the disorder is th en. , dp ne 
by a coherent potential approximation (CIPA).Ealla In 
previous work, Persson and Liebsch (PL)ll3 dealt with 
a similar model within the CPA. Their treatment ac- 
counted for non-retarded dipole-dipole interactions be? 
tween the particles. Later Stefanou and Modinos (SM)tj 
studied the same model. Their treatment of disorder 
effects used the simpler average T-matrix approxima- 
tion (ATA), however, they did include retardation effects 
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as well as sample-mediated interactions in their theory. 
Both PL and SM concentrated on calculating the aver- 
age transmission through, and reflection from arrays of 
relativelypSmall particles (i? ^ 10 nm). Meier, Wokaun, 
and Liaot£l have also studied dipole-dipole interactions 
in a self-consistent way in a disordered array of parti- 
cles by convoluting the response of a single-particle with 
a function describing the distribution of particles in the 

Our theory goes beyond earlier treatmentsliilli'0 in 
that it simultaneously includes effects of retardation, 
multipolar excitations, interactions with the sample, and 
disorder. Moreover, since we wish to calculate the in- 
tensity of diffusely scattered light, we do not only calcu- 
late disorder-averaged fields and particle polarizabilities 
(which one obtains from the CPA), but also disorder- 
averaged intensities. We calculate these intensities within 
a conserving approximation, thus using an averaging pro- 
cedure that includes exactly the same kind of scattering 
processes as the CPA calculation oi-the average fields. In 
the language of many-body theoryjHJ evaluating the aver- 
age polarizabilities is quite analogous to a single-particle 
self energy calculation, whereas the calculation of the dif- 
fuse scattering intensities requires that a vertex function 
(involving two-particle correlations) is solved for. 

The calculated results show good agreement with thg 
experimental ones. In particular, as in the experiment,El 
a series of resonance peaks emerge in the spectrum of dif- 
fusely scattered light at relatively small photon energies. 
A detailed study of the behind-lying mechanisms shows 
that these resonances occur whenever the field that drives 
the plasma oscillations in the silver particles reaches a 
maximum. This driving field is the sum of a contribution 
from the incident wave and a contribution from the wave 
refiected from the sample. The driving field displays a 
number of oscillations due to interference between waves 
refiected from the different interfaces in the multilayered 
sample. Thus, the interpretation we arrive at is simpler 
than the one proposed by SH.Ia It is, however, consistent 
with conclusions-drawn from earlier experimental results 
by Leitner et altn who studied light scattering from a sil- 
ver island film resting on a layered sample and observed 
characteristic changes in the scattering spectrum as the 
sample geometry was varied. 

The rest of the paper is organized in the following way. 
In Sec. |l| the basic theory is outlined, while Sec. Ill de- 
scribes the theory involved in the disorder-averaging. In 
Sec. I^the numerical results are presented and discussed, 
and Sec. |^ gives a brief summary of the paper. Three ap- 
pendices contain information on more technical aspects 
of the calculations. 



II. BASIC THEORY 

The system we consider consists of an array of spheres 
placed on or above a semi-infinite sample that may be lay- 



ered. A schematic illustration of the system is shown in 
Fig. |l|. The optical properties of the sample and particle 
materials are taken into account through local dielectric 
functions tabulated by PalikHj Since we will deal with 
relatively large particles non-local corrections to the di- 
electric properties should be relatively small. 

To begin with we assume that the array of spheres is 
perfectly ordered, but later we will relax this assumption 
replacing the array with a lattice gas where a particu- 
lar site is occupied by a sphere with probability p and 
unoccupied with probability 1 — p. 



A. Kirchoff integrals 

The basic task before us is to calculate the electric 
and magnetic fields everywhere in space given an incident 
wave with electric field E*"^', and magnetic field B'^''*. We. 
will use the vector equivalents of the Kirchoff integraH 
to do this. These formulae read (a time-dependence e~*"* 
is implicitly assumed everywhere) 

E(r) = E^'^^tr) + J dS' [tkc{n' x B(r'))G(r, r') 

+ {n' X E(r')) X V'G + (ri' • E(r'))V'G] (2.1) 



and 
B(r) 



dS' 



ik 



in' X E(r'))G(r,r') 



^{h' X B(r')) X V'G + [h' ■ B(r'))V'G] . (2.2) 



Here G denotes the Green's function of the scalar 
Helmholtz equation in free space, thus G can be writ- 
ten 



G(r,r') 



gifc|r-r'| 
47r|r - r'l 



and solves 

[V2 + A:2]G(r,r') --(^(^^(r-r'). 



(2.3) 



(2.4) 



wher e k — lu/c. The integrals appearing in Eqs. (2.1) and 
( ^.2|) run over all surfaces that enclose the nanoparticles 
as well as the sample, and the E and B fields appearing 
inside the i nteg rals ar e th e exact fields at these surfaces. 
Thus, Eqs. (2.1) and (2.2) are coupled integral equations 
from which the electromagnetic field always, at least in 
principle, can be calculated. 

The Kirchoff integrals become very useful also for de- 
tailed calculations once we have analytic expressions for 
the fields inside the spheres in the array and in the sam- 
ple. In the rest of this section, however, they will mainly 
be used as a tool to aid our thinking. 
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B. Multipole fields and sphere response 

The electromagnetic field inside and around each 
sphere can be expressed in terms of electric (E) and mag- 
netic (M) multipole fields. An electric pultipole field is 
defined by (we use Jackson's definitionsEEl) 



B(^) = -z,(fc,r)X, 



E 



-V X [zi{krr)Xim] , 



where kr is the magnitude of the wave vector kr — 
y/^^/c in a material with relative dielectric function 
Er-Ej A magnetic multipole is associated with the fields 

E(^^) = kzi{krr)X,rn; B^*^' = --V X [zi{krr)Xi,n] ■ 

c 

(2.6) 

In both these equations the vector spherical harmonic 



= {LYi„,)/ ^1(1 + 1). 



(2.7) 



is defined in terms of the angular momentum opera- 
tor L = —ir X V and the usual spherical harmonics 
Yim- The vector spherical harmonics fulfill the orthog- 
onality relations J X^;^/ (f2) • X;m(ri) = Sii>6mm', and 
/ dfiX;*,^, (fi) • [f X Xim(rj)] = on the unit sphere. 
The function zi{krr) stands for a linear combination of 
a spherical Bessel function ji^k^r), which is regular at 
r = 0, and a spherical Hankel function of the first kind 
hi{krr) = ji{krr) + ini{krr), which describes outgoing 
waves. 

In the present case, we can write the electric field inside 
the sphere at site i as 



E 



^kc[^\ii{krr)Xi, 



(2.8) 



Since the fields must be defined at the sphere center, only 
Bessel functions ji appear. Just outside the sphere, we 
have instead 



E = i: fc Ji{kr) + bl^PlHikr)] X,„, + 

{[4Z^Mkr)+b[Z^hi{kr)]X,„.]. (2.9) 



Irn 



-iW X 



Demanding that the usual boundary conditions for the 
E and B fields are satisfied at the surface of the sphere 
one can derive equations from which the field inside a 
sphere (the c coefficient) and the scattered field {b coef- 
ficient) can be calculated for each multipole, given the 
incident field (a coefficient). Since we focus on scattering 
properties, the key quantities are the "sphere-response 
functions" 



ekRj[{kR)+ji{kR) {e-l-Ji) 
e kRh[{kR) + hi{kR) [e-l-JiY 



and 



lAM) 

' AM) 
"■hn 



kR][{kR)-]i{kR)Ji 
kRh[{kR) - hi{kR)Ji' 



(2.10) 



(2.11) 



(2.5) where Ji is shorthand for 



Jl = krRj'i(krR)/jl{krR). 

Thanks to the symmetry of the sphere these response 
functions are independent of m. 

In the following, we will often describe the "state" of 
a sphere through the a and b coefficients. It is therefore 
convenient to collect these coefficients into vectors di and 
bi in "multipole space" with the structure 



(E) {E) 



^E) 



AM) 



AM) 



(2.12) 



etc 



,0 If 



we form a diagonal tensor V, in which the re- 



sponse functions sj^^'' and sj^""" enter in the appropriate 
places, the relation between incoming and outgoing waves 
at site i can be written 



AM) 



(2.13) 



C. Surface response 

Light will be scattered also from the sample surface 
and we need to deal with the surface response. A plane 
wave incident on the sample surface has a wave vector 



q- = qii - 2 



fc2 



(2.14) 



(the subscript "-" indicates that the wave propagates in 
the negative z direction) and can be written 



E 



2q_ -r 



(2.15) 



where the unit vectors for s and p polarization are s — 
[z X gy] and p = — [(?_ x (z x gy)], respectively. Here (j_ is 
a unit vector in the direction of the wave vector, and q\\ is 
a unit vector in the direction of the in-plane component 
of the wave vector. 

Th e reflected wave in response to the incident one 
( ^.15 ) can be written 



E 



refi - [xs (qil , oj)E'^'^s + Xp(q|| , Lo)E^^^p] e 



iq+-r 



(2.16) 
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where now q+ = qy + z^k"^ — \(\\\ p and p = — [g+ x (z x 
gy)] for this wave vector. In case the sample is homoge- 
neous the surface response functions Xs and Xv are found 
from the Fresncl formulae yielding 



(2.17) 



where ka is the wave number in the sample depending on 
the dielectric function of the sample es(w), 

kl = e,(u:)e. 

For p polarized light one gets 



fc2 ^ q-l - Jkl 



(2.18) 



Note that the expressions for all the surface response 
functions can be extended to the case of evanescent waves 
for which |q||| > k (when evaluating square roots the 
branch cut lies below the positive real axis). For lack of a 
better terminology, "plane wave" will sometimes be used 
to describe both propagating plane waves and waves that 
propagate in the directions parallel to the sample surface 
but are evanescent in the third, z, direction. 

When the sample consists of several layers of different 
materials the response functions Xv and Xs must be cal- 
culated by a generalized approach, for example a transfer 
matrix formalism, see Ref. Then one first makes an 
Ansatz for the electromagnetic field in the bottom layer 
of the sample in terms of one plane wave that propa- 
gates or decays exponentially in the downward direction. 
Then this wave is matched to a downgoing and an up- 
going wave at the bottom interface. These waves are 
propagated through the next layer, matched at the next 
interface, etc., until one reaches the sample surface where 
the ratio between the incident and upgoing wave ampli- 
tudes yields the surface response for the two different 
polarization types. 



D. Multiple-scattering solution for an ordered array 

Having introduced the necessary basic ingredients of 
the calculation we consider now the complete array of 
spheres (without the sample for the m oment), and let 
a plane wave like the one in Eq. (2.15) impinge on this 
system. Now Eq. ( [2.13| ) can be written 



(2.19) 



Thus the incident field at sphere i is, as can be under- 



stood by looking at Eqs . (2.1) and (2.2), a sum of one con- 
tribution, o^'"", coming from the projection of the plane 



wave in Eq. ( ^.15 ) onto the various multipoles, and con- 
tributions from waves scattered off all the other spheres. 
The explicit, rather lengthy, expressions for af"' and the 
coupling coefficients t^J (or rather its Fourier transform) 
are derived in Appendices ^ and ^ respectively. 

Thanks to the periodicity of the sphere array, t ij only 
depends on the re lative vector R — separating spheres 
i and j, and Eq. ( 2. IE ) can be solved by a Fourier trans- 
formation. For an arbitrary quantity Z defined on the 
lattice we introduce the Fourier transform 



Q 



and its inverse 



Q 



(2.20) 



(2.21) 



where N is the number of lattice sites. Applied to Eq. 
(^.19|) this yields 



5Q 



+ * Q Oq 



SO that 



5q 



(2.22) 



(2.23) 



When the sample is added to the problem one more 
surface, the plane z = zq — —R, will contribute to the 
Kirchoff integrals. The waves that are sent out from the 
sample surface are in turn emitted in response to the 
incident plane wav e or waves coming from the spheres. 
Consequently Eq. ( ^.13 ) now reads 



rref 



, (2.24) 



where "^''^ 



a^"' is the contribution from the incident wave 
after it has been reflected once from the sample, w is 
the sample-mediated self-interaction between the various 
multipoles on a sphere, while t is the sample-mediated 
interaction between two different spheres. These quanti- 
ties are discussed in Appendix |^. It should be kept in 
mind that both w and t describe events in which the 
source sphere sends out a wave that is scattered once off 
the sample and then goes directly to the receivi ng sp here; 
multiple scattering events enter through Eq. ( ^.24 ). We 
introduce a total driving field 



(2.25) 



and a total sphere-sphere interaction 

^ ^dir I ^sub 

ij ^ ij ' ij ' 



(2.26) 



and Fourier-transforming we find the solution for 6q 
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5Q 



1 — s w — s t 



Q 



so??*. 



It is also possible to rewrite Eq. (2.24) as 



->cxt 

a, 



t ij bj 



(2.27) 



(2.28) 



Thus, the sample-mediated self-interaction does not ap- 
pear explicitly in this equation, but is instead accounted 
for by the modified response tensor (1 — 's'w)~^*s. Inter- 
particle interactions mediated by the sample are of course 
still included in tij. The reasoning behind Eq. (2.28) in 
many ways makes it easier to think of the problem, since 
the sample in a sense has been eliminated from the for- 
malism. The effects of the sample on the spheres are sum- 
marized by the modification of the response tensor and 
the additional contribution if"*' to the particle-particle 
interaction. 



III. DISORDERED ARRAY 
A. Coherent potential approximation 

In the previous section we dealt with a system where 
each site on the lattice was occupied by a sphere. We 
will now treat a disorder model, namely a lattice gas. 
We still have a square lattice but only a random fraction 
p of the sites are occupied by spheres, whereas the rest 
of the sites are empty. 

The averaging over different disorder realizations will 
be done within the coherent-potential approximation.E3 
The CPA is known to be the best single-site mean-field 
theory, and it gives correct results in a number of im- 
portant limits. A certain realization of the disordered 
system is characterized by the response function of the 
sphere occupying the various sites; at the empty sites this 
response (V) vanishes. Within the CPA calculation we 
aim at finding an effective medium in which every site is 
occupied by average objects that can also be described 
by a surface response tensor that we will denote a in 
the following. As we saw above, the treatment can be 
simplified by incorporating the sample-mediated self in- 
teraction into the particle response. For the average par- 
ticles we do this by also introducing the tensor f3 which 
is related to a by 



/3 = ( 1 — aw 



<^ a={l+l3w) (i. (3.1) 



It should already now be pointed out that the ten- 
sors a and /?, describing the average scatterers, are not 
diagonal; the spherical symmetry is lost when interac- 
tions with the environment are taken into account. For 
example, in the case that the particles are treated as 
dipoles, the average particles have different polarizabili- 
ties in the directions parallel to, and normal to the plane 



of the array, respectively.E3 Both a and (i are, however, 
still lattice symmetric. Thus, only those tensor elements 
that couple multipoles belonging to the same irreducible 
representation of the array point group (in the present 
case where the particles interact with a sample, d^y) are 
nonzero. 

The response functions a and (3 are determined self- 
consistently by placing a real site, i.e. one which is occu- 
pied by a real sphere with probability p and empty with 
probability 1 — p in the effective medium of average par- 
ticles, and demanding that the real site does not cause 
any additional scattering on the average. 

To see what this means in detail, we consider the scat- 
tering off the real site as compared with when it is occu- 
pied by an average scatterer. An occupied real site yields 
the extra outgoing waves 



&1, extra — (/5oc " P) 



in a single scattering event. In Eq. (3.2) 



/3c 



-1 ^ 

1 — ¥w ) V, 



(3.2) 



(3.3) 



the effective response function for a real sphere includ- 
ing the sample-mediated self interaction, and ag could 
describe any combination of waves incident on the real 
site. In the CPA one includes not only single scattering 
events off the real site, but takes into account multiple 
scattering to all orders. This means that the extra scat- 
tering from an occupied site can be written with the aid 
of a T matrix Tqc, 

fcoc.oxtra = Toe OO = 1 - (/3oc - {(3 oc - P)a,o. 

(3.4) 

Here 700 is a propagator for multipole excitations de- 
scribing how waves sent out from the real site propagate 
through the effective medium of average scatterers and 
eventually return to the real, central site. One can show 
that 



700 = iV-i^(T-TQ/3 
Q 



t 



Q- 



(3.5) 



If instead the real site is empty, the extra scattering can 
be written 



^un, extra — ^un ^0 



1 - (Pun - P)lOO [Pun " l3)aQ. 

(3.6) 



Of course in this case the empty site does not give any 
scattering, i.e. = 0, and 



feun, extra = 1 + PjQO {- f3)ao 



(3.7) 
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The CPA self-consistency condition demanding that 
there should be no extra scattering on the average means 
that 



P &OC, extra + (1 - P) 6un,oxtra = 0. 



After some algebra Eq. (3.8) can be rewritten as 



1 



700 



(3.8) 



(3.9) 



Equation ( |3.9| ) has to be solved for /3 by an iterative pro- 
cess in which also 700 in Eq. (B.f;) has to be updated in 
each step of the iteration. 

From Eq. (3^) it is easy to verify one of the celebrated 
features of the CPA, namely that it approaches the cor- 
rect limit for both p = when (3 = 0, and p = 1 for which 
the self-consistent solution is /3 = /3oc = (1 ^ 'sw)^^'s. 



B. Diffuse scattering 

So far we have found the response functions a and f3 
that characterize an average effective medium that has 
the same effect as the lattice gas of particles on the in- 
cident light in terms of coherent transmission and reflec- 
tion. 

To calculate the intensity of diffusely scattered light 
into a certain direction, the theory has to be developed 
further. Up to now we have only calculated average quan- 
tities involving one single multipole excitation, but to 
find diffuse scattering intensities averages involving two 
multipole excitations must be calculated. 

Suppose that a plane wave described by the ampli- 
tude vector flq' impinges on the disordered particle ar- 
ray. In response to this there will of course be multi- 
pole excitations at the incident wave vector Q that can 
be characterized by 6q. But as a result of the disorder 
there are also multipole excitations described by 6q' at 
other wave vectors Q'. Now, averaging over a large num- 
ber of different disorder realizations (in our case, lattice 
gas arrangements) would show that (6q) 7^ 0, (here {X) 
denotes disorder- averaging of the quantity X) whereas 

(6q') = 0, since the diffusely scattered light has an es- 
sentially random phase that varies from one disorder re- 
alization to another. Yet averages like (6J q'^C.Q') ^^'^ 
C are shorthand indices for different multipoles) of course 
do not vanish and hence there is diffusely scattered light. 

Assume that we have managed to calculate the coeffi- 
cients b(^Qi for one particular realization of the disorder. 
The (reflected) scattered electric field Esc (r) , with a wave 
vector 



fc2 



q'll +q'J with q'^ 



whose in-plane wave vector q'|| = Q'-l-G equals Q' up to 
a reciprocal lattice vector G = {2t: / a){nxX -\- nyy) with 
rix and riy being integers, can then be written 



Esc(r) 



(2^ 



Here the function 



rf?(q') =5^(q'+) + e"^V^I^°l;^,(q'll).g^(q'_) (3.11) 

describes coupling from multipole C to an outgoing s po- 
larized wave, either directly (first term) or through a re- 
flection off the sample (second term). The functions 5^'^ 
which are defined in Appendix ^ measure the "overlap" 
between multipole C and a plane wave with wave vector 
q'_i_ = q'll ± q'^z. When the wave leaves the system with- 
out being reflected off the sample the wave vector should 
have a positive z component, whereas when a reflection 
takes place the multipole must first couple to a down- 
going plane wave (q'_), and the amplitude of the wave 
is also affected by a phase shift due to the additional dis- 
tance that the wave must travel and the surface response 
function Xa- In the same way. 



d^(q') = .<(q'+) + e'^V^I^°l;,^(q'||)g^(q'_) (3.12) 

for outgoing p polarized waves. The scattered power is 
now for this particular realization 



(Pq' 



d«(q')*ci^(q') + rff(q')*<(q') 



(3.13) 



If we were to calculate the scattered power for another 
realization of the disorder we would arrive at the same 



formal expression (3.13); the functions <f^^\c() would 
still be the same, however, the values of the h coefficients 
would of course be different. Thu s, the disorder-averaged 
scattered power is given by Eq. ( 3.13 ) provided that the 
product 6|q/^cQ' replaced by its disorder-average 

Obviously, to find the intensity of diffusely scattered light 
we must calculate the averages (&Jq/&cQ')- 



C. Vertex function 

Since we have calculated the properties (essentially /3) 
of the average objects within the CPA, also the calcula- 
tion of (6|q/ ^^CQ' ) niust be done at the same level of ap- 
proximation in order to conserve the total energy. This 
means that this calculation must take into account re- 
peated single-site scattering events to all orders. To this 
end we have to solve a Bethe-Salpeter equation for a ver- 
tex function A tjlia,t sums all "ladder diagrams" as illus- 
trated in Fig. 
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The irreducible vertex S describes the simplest scatter- 
ing event, off one real site, that enters the sum of ladder 
diagrams, ft is given by an average between the scatter- 
ing intensity from an occupied site and from an unoccu- 
pied site, and can be expressed in ter ms of t he p reviously 
introduced T matrices [see Eqs. (3^) and (|3.6| )] as 



Sfi =p(r:Jc,(Toc)c5 + {l-p){T*MTu^)c5. (3.14) 

The full vertex function A is then found by allowing 
for the possibility of consecutive scattering events off dif- 
ferent sites. Summed to all orders this yields the Bethe- 
Salpeter equation (summation over repeated indices is 
assumed implicitly) 



A 4'; — I rr'''''" a ''"'' 



M"5 ■ 



(3.15) 



In Eq. (3.15), 11 describes the propagation of the multi- 
pole excitations through the effective medium from one 
scattering site to the next. This propagator involves es- 
sentially the same quantities as did the "return" propa- 
gator 7oo encountered earlier. We can write 



n 



fs = "(7oo)c^(7oo)c5 + }_^(7Q)e.,(7Q)c5, (3-16) 



1 V^i- 



where 



7Q 



1 



Q 



-1 

tQl3] t 



Q- 



(3.17) 



I nsert ing the above expressions for S and 11 into Eq. 
( |3.15 ) we can solve for A numerically by means of a ma- 
trix inversion. 

Finally, we must relate the vertex function to the av- 
erages (fe|q/6cQ') needed to find the diffuse scattering 
intensity. One can write 



Here 



xA«;^ (a<=«*),((f«),. 



Q 



(3.18) 



(3.19) 



represents the effective driving field at the site where the 
first scattering event occurs. This effective field results 
from screening of the external field (given by Oq*) by the 
average particles in the effective medium. In the same 
way, th e ten sor (1 — /3iQ')^^ appearing in front of A 
in Eq. ( 3.18 ) accounts for the screening of the outgoing 
waves by the average particles after the last scattering 
event has occurred. One should note that the outgoing 
wave vector Q' occurs in this expression, whereas, in the 
case of an incident plane wave, there is only one term (at 
wave vector Q) in the sum defining a^^ . 



IV. RESULTS AND DISCUSSION 

Let us consider a system of spheres with radius i? = 50 
nm and lattice parameter a = 120 nm. These values 
probably describe the average properties of the particle 
overlayers studied by Stuart and HalH rather well. How- 
ever, one should keep in mind that in that experiment 
the typical particle shape was most likely not spheri- 
cal. Judging from the manufacturing process, the par- 
ticles studied in that experiment looked more like oblate 
spheroids than spheres. 

Figure || shows calculated results for the intensity of 
scattered light from a disordered array of silver particles 
placed on a Si/Si02/Si sample. In close agreement with 
the experimental results found in Ref. || the scattering 
intensity in Fie. ^ (a) has two isolated peaks at 1180 
and 710 nm.EB These isolated peaks are then followed 
by a broader "scattering band" at shorter wavelengths 
^ 550 nm. We also note that the overall intensity of the 
scattered light is in reasonable agreement with the exper- 
imental results; the measured diffuse scattering intensity 
collected over a solid angle of « 27r/5 was about 1 % of 
the incident intensity, while we here get intensities of the 
order of 10 % when calculating the light collected over a 
solid angle of 27r. 

What mechanism makes the two isolated, long- 
wavelength peaks appear in the spectrum? These fea- 
tures only appear when there is a layered sample; as can 
also be seen in Fig. ^ (a) there are no such peaks in 
the spectrum calculated with a homogeneous dielectric 
sample (LiF). Stuart and HallQ argued that the inten- 
sity of the scattered light reaches a maximum whenever 
the particle-particle interaction mediated by waveguide 
modes in the sample reaches a maximum. While this 
certainly is an interesting suggestions the results of our 
calculation show that the isolated resonances have a sim- 
pler explanation. 

The isolated peaks in the scattering intensity coincide 
with maxima in the electric field that polarizes the metal- 
lic particles in the first place. The dotted curve in Fig. 
^(a) shows the squared strength of the driving electric 
field found in the plane of the sphere centers as a re- 
sult of the incident wave and its reflection off the sample. 
This field has been calculated in the absence of the parti- 
cle overlayer. In the present situation, when the incident 
field hits the sample at normal incidence (no physical dif- 
ference between s and p polarization), the driving field 
can be written 

Edrivc = i?^^^s[l + e2^'=l^"lxs(q|| =0,c.)], (4.1) 

thus the quantity plotted together with the scattering 
intensity in Fig. is 



1 



Xs(q|| = o,w) 



The oscillations in the driving field strength with varying 
wavelength is consequently a result of constructive and 
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destructive interference between waves reflected from dif- 
ferent interfaces of the multilayer sample. For the case 
shown in Fig. || (a), even if the presence of the other lay- 
ers also affects the driving field, the thickness di, of the 
top Si layer is the most important parameter in determin- 
ing the interference maxima and minima; Agi = X/y/es[ 
equals 2di when Tiui « 1.1 eV (A wll30 nm), and 4di 
when huj « 1.9 eV (A «650 nm). Of course, when one 
changes the sample geometry i.e. the thicknesses or the 
composition of the various layers, the maxima of the driv- 
ing field shift in frequency. As can be seen in Fig. ^ (b), 
where the thickness di of the top Si layer was increased to 
240 nm, this brings with it identical shifts of the isolated 
peaks in the scattered intensity. This further supports 
the interpretation given here. 

Figure ^ shows how the in-plane electric field varies 
with z in the Si/SiO/Si sample used for the calculations 
in Fig. 1^ (a) when there is no particle overlayer and waves 
with photon energies 1.1 eV, 1.3 eV, 1.5 eV, and 1.7 eV, 
respectively, impinge on the surface at normal incidence. 
As can be seen there are reflection resonances when it is 
possible to fit in a wave that has antinodes at the inter- 
faces of the top Si layer. 

A simple explanation for why surface waveguide modes 
are not causing the scattering resonances lies in the very 
small scattering cross section for a single sphere at the 
relevant wavelengths. This cross section is of the order 
{kny R^. With a wavelength of 1000 nm and i? = 50 nm 
kR ~ 0.3, which means that the scattering cross section 
is of the order of 1 % of the geometric cross section. 
Now, the process discussed by Stuart and HallO involves 
two scattering events off nanoparticles because the in- 
coming light must first be scattered into a waveguide 
mode by the particles and then out again in order to 
be observed as diffusely scattered light. Of course such a 
process will not be very effective, unless some other fac- 
tor compensates for the small scattering cross sections. 
In this case the waveguide modes do offer some compen- 
sation; the response of the semiconductor surface is very 
strongly enhanced for combinations of in-plane wave vec- 
tor qy and photon energy hoj that coincide with those of 
a waveguide mode. But this only happens for a narrow 
interval of wave vectors. There is not enough phase space 
for the waveguide-mode-mediated particle-particle inter- 
action to become an important factor for the scattering 
of long- wavelength light. 

This point is best illustrated by looking at a plot of 
the particle-particle interaction as a function of in-plane 
wave vector shown in Fig. ^. The curves show the value 
of |/3ii,iiQiii,iiQ| as a function of Q for two different 
photon energies. This quantity yields the ratio between 
the secondary wave hitting a sphere as a result of a pre- 
vious scattering event off the other average particles, and 
the primary incident wave. With this particular choice 
of tensor elements we are focusing the attention to the 
waves sent out by in-plane electric dipoles on the par- 
ticles. We note that for most wave vectors the sphere- 
sphere coupling for these photon energies lie well below 1, 



i.e. the primary wave dominates in strength over rescat- 
tered secondary waves. For certain wave vectors the 
sphere-sphere coupling is, however, rather strong. This 
happens when the electromagnetic oscillations of a parti- 
cle are in resonance with a waveguide mode of the sample. 
In this case the sample-mediated coupling has a sharp 
peak. With huj = 1.0 eV there are two such resonances. 
The peak marked A' is due to an s polarized waveguide 
mode whereas the peak B' is caused by a p polarized 
mode. When increasing the photon energy these peaks 
move out towards higher Q. At 1.3 eV their correspond- 
ing positions are given by the peaks A and B. The new 
peak C is associated with a standing wave resonance very 
similar to the one occurring at Q = at 1 eV. 

In this context we should also say that the waveguide 
modes may play a much more pronounced role in exper- 
iments with periodically ordered arrays of nanoparticles. 
In this case the surface modes that can be excited must 
have an in-plane wave vector that either equals the inci- 
dent one or else differs from it by a reciprocal wave vector 
G. In that case the sample-mediated particle-particle in- 
teraction is not averaged over different wave vectors in 
the same way as in the present case. The recent experi- 
ment by Linden et aln who measured light transmission 
through a layered sample covered by an array of gold 
nanoparticles, provides such an example. 

Turning to the range of the spectrum at higher photon 
energies, the features appearing there are caused by a 
combination of particle-sample and particle-particle in- 
teractions. The disorder in the sphere array of course 
leads to considerable broadening of the spectrum. To 
analyze this in more detail, we display in Fig. ^ (using a 
photon-energy scale) the diffuse scattering intensity for 
three different coverages p =0.1, p =0.2, and p =0.5, re- 
spectively. It is evident that the particles do not act in- 
dependently of each other already from the fact that the 
scattering intensity as a function of coverage saturates 
for photon energies above 2 eV or so. Overall, p = 0.2 
gives a stronger scattering than p — 0.5. This is in con- 
trast with the situation at the peaks at 1.1 eV and 1.7 eV 
where the scattering intensity is basically proportional to 
p, consistent with a picture where the different particles 
scatter light essentially independently of each other. 

While the effects of particle-particle interactions are 
less clean-cut in the previous case than when the parti- 
cles are much smaller than A/27r (cf. Ref. |l^), we can still 
gain some more insight into what is happening in the var- 
ious parts of the "scattering band" by looking at the field 
distribution axpund the average objects at different pho- 
ton energies.tZl Figure [t] shows the electric field around 
an average scatterer at 2.5 eV, 2.9 eV and 3.6 eV, when 
an array of such scatterers are placed on a Si/Si02/Si 
sample with di = 160 nm and d2 = 200 nm and illu- 
minated by light at normal incidence. The solid arrows 
represent the real part of the field (in phase with the 
incident wave in the plane z = 0). (i) At 2.5 eV there 
is a concentration of the field to the region in between 
the particles. This is a result of direct particle-particle 
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interactions where a particle is polarized by the effective 
field which is the sum of the external field and the field 
from the other particles, the nearest neighbors giving the 
largest contributions. Indeed, if one considers two (real) 
silver spheres with R — bO nm at a distance 120 nm from 
each other that are driven by an external field directed 
along the line joining the two spheres, they have a col- 
lective resonance at about 2.6 eV. Thus this resonance 
is redshifted relative to the dipole resonance of an iso- 
lated particle which occurs at about 3.1 eV if i? = 50 nm 
since the field from the neighouring sphere acts in phase 
with the external field when the frequency lies below the 
single-sphere resonance frequency. In Fig. ^ we see that 
as the coverage is increased, the low-frequency edge of 
the scattering band shifts to the left and the intensity in 
this part of the spectrum increases, consistent with the 
growing importance of particle-particle interactions at a 
larger p. (ii) Looking at the fields at 2.9 eV we see that 
they are more concentrated to the region in space near 
the sample. Here particle-sample interactions play a rela- 
tively important role. What happens is that the external 
field polarizes the particles, and the particle dipole mo- 
ments, in turn, induce image charges in the sample that 
reinforce the driving electric field. In this way the dipole, 
but also higher multipole moments on the particle are fur- 
ther amplified, leading to the field concentration near the 
sample surface, (iii) At 3.6 eV, on the other hand, the 
fields are concentrated to the upper (vacuum) side of the 
particles, and the fields are also weaker in this case than 
in the previous ones. Here the induced dipole moment on 
a sphere is essentially out of phase with the driving field. 
Hence the image charges in the sample this time set up 
an induced field opposing the external electric field acting 
on the particle. The enhancement mechanisms discussed 
above are no longer efficient and the fields are localized 
to the top side of the particle. j-. 

In their experiments, Stuart and Halld also studied 
light scattering from an overlayer of Ag particles on a 
silver sample covered by spacer layers of LiF of varying 
thickness. The silver surface possesses elementary exci- 
tations in the form of surface plasmons and these exci- 
tations can mediate interactions between the silver par- 
ticles in the overlayer. Also in this case, the experimen- 
tal results exhibited scattering resonances, albeit very 
broad ones, at rela|t|ively small photon energies. In the 
experimental papeitl they were discussed along the same 
lines as the resonances ocurring in the case of the layered 
Si/Si02/Si sample, namely as the result of strong surface- 
mediated particle-particle interactions that enhance the 
scattering intensities. 

Figure g shows results for the calculated scattering in- 
tensity from a particle overlayer on silver samples with 
LiF spacer layers of two different thicknesses, 120 nm and 
200 nm. With a Ag sample the calculations are numer- 
ically more difficult than with a Si sample because the 
surface response now, due to surface plasmons, has con- 
siderably sharper features than what one sees in Fig. ||. 
There is therefore some numerical noise in the results, as 



can be seen as irregular oscillations of the curves in Fig. 
I- 

The results in Fig. ^ agree well with the experimental 
ones from a qualitative point of view. The shape of the 
spectra are the same as in the experiment, even if the 
positions of the maxima are somewhat shifted compared 
with the results in Fig. 4 of Ref. ^ Part of the expla- 
nation for this may be that the height of the particles 
is smaller in the experiment than in the calculation. As 
was the case with a Si/Si02/Si sample, the calculated in- 
tensity of the scattered light closely follows the strength 
of the driving field, which is also plotted in the figure. 
We can conclude that also for the silver sample the over- 
all shape of the diffuse scattering spectrum is set by the 
driving electric field which in turn is affected by inter- 
ference effects between waves reflected from the different 
interfaces in the sample. Increasing the thickness of the 
LiF layer from 120 nm to 200 nm shifts both the peaks 
in the scattering spectrum towards longer wavelengths. 
This brings with it a decrease in the scattering inten- 
sity at long wavelengths since the maximum of the driv- 
ing field falls at longer wavelengths where the scattering 
cross section of the particles is smaller. The scattering in- 
tensity at the short-wavelength peak, on the other hand, 
increases when d increases from 120 nm to 200 nm. Both 
these trends for the scattering intensity are qualitatively 
consistent with the experimental results. The strength 
of the short-wavelength peak grows when d changes from 
120 nm to 200 nm because the maximum of the driving 
field in the latter case coincides with the previously dis- 
cussed resonance near 2.5 eV at which the in-plane dipole 
moments of neighbouring particles oscillate in phase [cf. 
Fig. 0(a)]. 



V. SUMMARY 



In this paper we have presented a formalism by which 
the diffuse light scattering intensity from a disordered 
array of spherical nanoparticles can be calculated. The 
theory is built on an expansion of the electromagnetic 
field around the nanoparticles in terms of multipoles, and 
takes into account interparticle interactions mediated by 
near and far fields as well as the sample. The theory was 
then applied to systems similar to those studied exper- 
imentally by Stuart and HalB in which silver nanopar- 
ticles are placed on layered substrates that support var- 
ious surface waves. The results of the present calcula- 
tion agrees rather well with the experimental ones. In 
particular the intensity of diffusely scattered light from 
nanoparticles on a layered Si/Si02/Si sample shows sev- 
eral characteristic peaks at long wavelengths. We have 
shown that these peaks are due to reflection resonances 
of the sample and not enhanced particle-particle interac- 
tions mediated by sample waveguide modes. 
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APPENDIX A: COUPLING TO PLANE WAVES 
1. Plane- wave-to-multipole coupling 

In this Appendix we present results for the coupling 
factors between multipole fields and plane waves. For a 
detailed derivation we refer to Ref. 

We consider first a plane wave written on the general 
form, 



E 



(r) = ^^E^'^s + E^P^pj e^'i"" 



(Al) 



where s = [z x q||] and p = — [q x (z x q||)] that is inci- 
dent on a sphere. This wave can be expanded in terms 
of electric and magnetic multipoles around the center Rj 
of sphere j as 



E(r) = J2k jiikr)Xi„, + zV x U^\n{kr)Xirn 



Im 



(A2) 



where 



|r - R, 



The coefficients in this expression depend linearly on 
the amplitudes of the incoming plane wave, i.e. 



(E) _ fEp 



f^ME^^^ + fE:{ci)E^'\ and 



where 



and 



The quantities, ?7;m(q) and V;„i(q), are given by 



(A3) 
(A4) 
(A5) 



+s_F_(?,TO)l^,_„,+i(r!g)], 



(A6) 



where Qq denotes the direction, the angles 9q and (/)q, of 
q and 



s± — s ■ X zL is ■ y, 



and 



+T]^F_{l,m)Yi^^^+i{nq) -?7,2my,,_™(r!,)], (A7) 



with 



and 



fi = k ^qx(iX(7||); r]± = ij^ ± iijy , 



Fj^{l,m) ~ \/(/ — m)(Z + 771+1), and 

F_(/, 777) = ^{l+7n){l — 777 + 1). 



(A8) 



The above expressions are also valid for evanescent 
waves that decay exponentially in the positive or negative 
z direction. The most general wave vector we consider 
here lies on the mass shell, and can be written 



q = q|| ± z ./P _ |q 



(A9) 



The plus sign holds for waves that are outgoing (propa- 
gating or evanescent) in the positive z direction, and then 
cos 6q is either real and positive, or lies on the positive 
imaginary axis. The minus sign holds for outgoing waves 
in the negative z direction. 



2. Multipole-to-plane-wave coupling 

We also need to know how much the outgoing waves 
from a multipole contributes to the amplitude of a certain 
plane wave. To this end we apply the Kirchoff integrals. 
Thus, for a magnetic multipole we calculate E from Eq. 
(2.1), while for an electric multipole we calculate B from 
Eq. (^1^) . To carry out the integration we use the follow- 
ing form for the Green's function 



G(r,r')=7 



.iqir(r||-r ii) 



(AlO) 



The scalar plane wave appearing here is expanded in 
terms of spherical harmonics, so that the overlap inte- 
grals with the multipole fields can be carried out. As 
is evident from Eq. (AlC), the radiated field is a linear 



combination of many plane waves, and the result can be 
written as 



E(r)= J^a^l 



frm(q)p]e 



zq-r 



(All) 



The coupling factors are found to be 



5L!(q) = 5/t'(q) 



2ypHqi 



(-l)'+'"+'C/,,_„(q), 
(A12) 
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and 



(-1) 



2yF 



-v^,-m(q)- 



(A13) 



APPENDIX B: CALCULATION OF 
INTERPARTICLE COUPLING 

The Fourier transform of t is given by 



*7*dir _ Vdir -iQ (Ri 

t Q - Z^'^ij ^ 



(Bl) 



where i can be any site. If aU spheres in the array sends 
out waves described by 6e**^'^J , the waves received at site 
i (without any intermediate scattering taking place) is 

= e*'*-^' T^'' b. (B2) 
In particular, at the site at the origin, 



ao 



t dir ; 

Q ' 



We therefore calculate fQ by identifying the waves inci- 
dent on the sphere at the origin. 

We can simplify the detailed calculation by realizing 
that outside a source sphere the outgoing waves appears 
to originate from the center of the sphere, and likewise, 
when expanding the waves incident on the sphere at the 
origin in different multipoles, this can be done at any dis- 
tance from the sphere center. In the calculations at hand 
here we therefore set the radius of the source sphere to 
r2 and that of the receiving sphere to ri and let both 
these radii tend to zero. The physical size (the radius R) 
of the real spheres is of course relevant to the physics, 
but this enters the calculation only through the response 
functions s|^^ and Si^\ 

Consider now an electric multipole source with angu- 
lar momentum quantum numbers I' and m' . We calcu- 
late the B field this generates near the origin using the 
Kirchoff integral in Eq. (pT^), where now the last term 
vanishes, so that 



B(ri) = 5: 



ik 



^2 / ^^^2 

+ {f2 X B(r2))V2G(ri,r2 



(f2 xE(r2))G(ri,r2-Rj) 



(B3) 



Here the phase factor e**^'^^ has been included explicitly, 
and B(r2) and E(r2) should be evaluated on the sphere 
at the origin (even if this, due to the condition j ^ 0, 
does not contribute to the sum). 



The summation over lattice sites only involves the 
Green's function and the last phase factor, two scalar 
quantities. We therefore change the order between in- 
tegration and summation and begin by evaluating the 
lattice sum 

5Q=^G(ri,r2-R,)e*'^-^^ = 



y-v exp [ik\ri — r2 — R 

3^0 



47r|ri - r2 - R 



(B4) 



■J I 



and its gradient W2SQ. The calculation uses Ewald 
methods borrowed from KKR theoryH The sum can be 
expressed both in terms of a one-center expansion 



Sq 



00 L 

E E 

L=Q M=-L 



Dlm 3 L{kri2)YLM {^12), 



(B5) 



where ri2 and rii2 denote the length and direction of the 
vector ri2 = ri — r2 , and a two-center expansion 



•Sq = ^ Aijni'm.'3i{kri)ji>{kr2)Yim{^i)Yi*^,{^l2) 



lull' 



(B6) 



The two sets of expansion coefficients are related through 



A 



imi'm' = 47r i' ' ^ Dlm Cf"^ 

LM 



Irnl'ni' ' 



where Cf"J^l^, denotes a Gaunt integral 



G 



LM 
Iml'm' 



dnYi*„^{n)YLM{^)Yi,m'{n). 



(B7) 



(B8) 



The two expansions serve different purposes. When eval- 
uating B from Eq. (B3), and then projecting the result 
onto different multipoles at the origin, we will use the 
two-center expansion. The one-center expansion is use- 
ful because using Ewald methods Sq can be split into a 
long-range part and a short-range part, and the coeffi- 
cients Dlm can be calculated rapidly. 
There are three contributions to Dlm, 



Dlm^D^j^I + D^j^I + D^^Slo. 

1(1) 



(B9) 



The long-range contribution Dj^^^ can be expressed as a 
sum of crystal-rod integrals 



D 



(1) _ 

|q 



G 



X ■ 



G|2 



■ ^LM(^q-l-G) 



,(fe^-|q+G|^)/,, 



, (BIO) 



1(2) 



here q = Q + zq^ and 77 is a separation parameter. D)^^^ 
comes from a direct summation over the lattice, and thus 
collects short-range interactions 
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2L+1 



(Bll) 

where the integral 



(B12) 



can be related to the complementary error function. The 
last contribution compensates for the fact that D^qq con- 
tains some self-interaction contributions, 



D. 



(3) 
00 



ik 



27r^n!(2n-l) 

^ n— ^ ' 



(B13) 



In the calculations here the separation parameter 77 was 
related to the lattice parameter a as 77 = 3.24/a^; the 
value of T] of course does not affect the result for Dlm- 

To proceed with the calculation of tq, we insert the 
two-center expansion Eq. (B6) into Eq. (B3), together 
with expressions for f2 x E and f2 x B for a certain elec- 
tric multipole {L'M'). The resulting surface integrations 
are now straightforward since Xl' m' is a linear combina- 
tion ofYL,^M'+i, Yl'^m', and Yl'^m'-i- We get 

X {n-F+{L', M') Ai,^x'M'+i 

+h+F- {L', M') Ai^m,L',M'-i + n,2M' Ai^^x',M' } • (B14) 

In this expression n+ ~ x ± iy, fiz — z, and and F- 
are defined in Eq. ( A8). 

Next we project B onto different multipoles around the 
origin. The field corresponding to an electric multipole 
{LM) can be written 



k 



B = '-a^LMjL{kri)'^LM{^i)-, 



therefore 



7 n 

-afljL{kr^) = j df7i(XiM(f^i))* • B(ri), (B15) 



and the final result for t 



(EE) 

LML'M'Q 



'M'Q ^ '^lm/^l'M' reads 



(EE) 



1 



LML'M'Q 2fc ^L(L + 1)L'(L' + 1) 
X {F+(L, M) F+{L', M') Alm+i,L',m'+i 
+F_ (L, M) M') Ai,M-i,L',M'-i 

+2MM'Ai,M,L',M'}- (B16) 



The corresponding calculation for the projection onto 
a magnetic multipole (LM) is somewhat more involved. 
We use the fact that for a magnetic multipole 



1 



ri • B(ri) = -a[^l)jL{kn)^LiLTT)YLM{^i), 



and consequently 
-a^LMjL{kri)^L{L 



(B17) 



Here the Cartesian components of ri introduce spherical 
harmonics Yn, Yio, and Yi_i into the integral. The final 
result thus contains Gaunt integrals, 



{ME) 
LML'M'Q 



2L + 1 



2k ^L(L + l)L'{L' + 1) 

^LM,L~l,M+l ^+(^'^ ^'^') ^L-l,M+l.L' M' + l 



'^'^LM,L-1.M-1 , M') Al_i^m-1,L' ,M' -I 



1,0 

LM,L-1,M 



2M' Al^ 



1,M,L',M' 



} 



(B18) 



We have determined half of the elements in Iq. If the 
sources instead are magnetic multipoles we use the other 
Kirchoff integral in Eq. ( ^.l| ) to calculate E(ri) around 
the origin. The calculations are perfectly analogous to 
the ones above and in the end one finds the following 
symmetry of t , 



(MM) 
LML'M'Q 



t 



{EE) 

LML'M'Q^ 



and 



{EM) 

LML'M'Q 



{ME) 

LML'M'Q- 



(B19) 



APPENDIX C: CALCULATION OF 
SAMPLE-MEDIATED INTERACTIONS 

Here we will calculate the sample-mediated self- 
interaction w as well as . 

Starting with w, suppose a certain multipole {I'm! a') 
is excited on a sphere. Then plane waves are sent out 
from the sphere, reflected off the sample surface either as 
s or p polarized light, and after that the wave impinges 
on the sphere again. The strength of the incident 
wave can be calculated using the previously introduced 
functions / and and the result reads 



Iml' m' 



(2^ 



- /r™ " (q+)x." (qil )3;C'' (q_ V^l^°l , (CI) 



where 



q± = qii ± zJk"^ 



12 



and a" denotes the plane- wave polarization (s or p) . The 
angular integration (over the different directions of qy) 
is straightforward, and yields a finite result only when 
m = to'. The remaining integration over |q|[| must be 
done numerically. 

Once w is known, it is relatively easy to calculate t q^. 
If the same multipole is excited on all the spheres, with 
relative phase e'*^ -^^ , the incident wave at the origin 
must be proportional to a sum of t and w. This sum 
is easily calculated thanks to the periodicity of the radi- 
ating system. A certain tensor element is found through 
a summation over the 2D reciprocal lattice 



^ccll ^ 

(7 



iSub.CTCr' 

''I 



/r™ " (q+ )Xa" (qil )ftC'' (q- )e''^ V^l^«l , (C2) 



where now 



Q + G, and q± 



±zJk 
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Johansson, FIG. 1 
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FIG. 1. Schematic illustration of the sample and overlayer of spheres. The two sample layers have thicknesses d\ and 
di, respectively. A lattice (lattice gas) of spherical particles with radius K and lattice parameter a covers the sample. The 
particles interact through direct t fj"^ and sample-mediated t interactions. (The empty space between the particle layer and 
the sample is there only for illustrative purposes; the particles do rest on the sample.) 
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Johansson, FIG. 2 

FIG. 2. Diagrammatic representation of the Bethe-Salpeter equation illustrating how the vertex function corresponding to 
the diffuse scattering intensity is built up by contributions from a single scattering event (first term) as well as by contributions 
from repeated scattering events. 
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Johansson, FIG. 3(a) 
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FIG. 3. (a) Relative intensity of the diffusely reflected light (integrated over the full solid angle 27r) from a disordered 

ovcrlaycr of Ag spheres. The sample is illuminated at normal incidence. The two solid curves show results obtained with a 
layered Si/Si02/Si sample and with a homogeneous LiF sample, respectively. The dotted curve shows the magnitude (squared) 
of the driving electric field relative to the incident field in the plane of the sphere centers. This field was calculated in the 
presence of the Si/Si02/Si sample, but without the Ag particles, (b) The corresponding results obtained with an increased 
thickness of the top Si layer. The maxima in the driving field have shifted to larger A [the 1200-nm-peak from (a) falls outside 
the figure at w 1600 nm] and the scattering-intensity maxima follow. 



15 



Johansson, FIG. 4 
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FIG. 4. The magnitude of the in-plane electric field as a function of z in the absence of any particle overlayer when a wave 
with amplitude 1 impinges at normal incidence on a Si/Si02/Si sample with the same geometric parameters as in Fig. ^ (a). 
Results are shown for four different photon energies as indicated next to the curves. The corresponding free-space wavelengths 
are 1130 nm, 950 nm, 825 nm, and 730 nm, respectively. Note that the value of the field magnitude at z — 0, the sphere center, 
equals the driving field whose square is plotted in Fig. ^ (a) . 



Johansson, FIG. 5 



D) 0.75 



O 
O 



a. 




0.25 



r M 
In-plane wave vector 

FIG. 5. The sphere-sphere coupling, or more precisely the product l/3ii,iiQtii,iiq|, plotted for two different photon energies 
1 eV (thick curve) and 1.3 eV (thin curve) as a function of the in-plane wave vector Q along the diagonal of the Brillouin zone 
from the center point (F) where Q = to the corner point (M) where Q = {2n/a){x + y). The peaks marked A and A' are 
due to s polarized waveguide modes the ones marked B and B' are due to p polarized waveguide modes. The curves exhibit 
discontinuities (coming from the interaction t) at the point where |Q| = k. 
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Johansson, FIG. 6 
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FIG. 6. Diffusely reflected light intensity from a Si/Si02/Si sample when the coverage p of Ag spheres is varied as indicated 
next to the curves. The other parameter values are the same as in Fig. (a). 
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Johansson, FIG. 7(b) 
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FIG. 7. The fields around an average scatterer for three different photon energies: (a) 2.5 eV, (b) 2.9 eV, and (c) 3.6 eV. 
The parameter values are the same as in Fig. I (a). The full arrows give the real part of the field (in phase with the incident 
wave) and the broken arrows the imaginary part. The arrow in the middle of each figure shows the strength of the field of the 
incident wave [thus, the scaling used in panel (c) differs from that of the other panels]. 
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Johansson, FIG. 8(a) 
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FIG. 8. Relative intensity of the diffusely reflected light (integrated over the full solid angle 27r) from a disordered overlayer 

of Ag spheres with a relative coverage p = 0.5 on a LiF/Ag sample. In both plots the full curves show the scattering intensity, 
whereas the dotted curves show the result for the (squared) driving field. The thickness of the LiF layer is (a) 120 nm and (b) 
200 nm, respectively, in the two plots. 
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